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We consider integrable quantum spin chains with alternating spins 
(Si, 52). We derive a finite set of non- linear integral equations for the 
thermodynamics of these models by use of the quantum transfer ma- 
trix approach. Numerical solutions of the integral equations are pro- 
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vided for quantities like specific heat, magnetic susceptibility and in 
the case S± = S% for the thermal Drude weight. At low temperatures 
one class of models shows finite magnetization and the other class 
presents antiferromagnetic behaviour. The thermal Drude weight be- 
haves linearly on T at low temperatures and is proportional to the 
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central charge c of the system. Quite generally, we observe residual 
entropy for Si ^ S2. 

PACS numbers: 05.50+q, 02.30.IK, 05.70Jk 
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1 Introduction 



Integrable quantum systems and their associated classical vertex models have 
been extensively studied in the last decades [HE]- A large part of these 
systems is exactly solvable by Bethe ansatz techniques providing spectral 
data and in some cases also the eigenvectors. 

After establishing the integrability and deriving the exact solution for 
the spectrum, the main questions one likes to answer concern the physical 
properties of the system in dependence on temperature, magnetic field etc. 
There are many investigations of integrable system in the thermodynamical 
limit at finite temperature. In fact, we have several established routes to 
this goal. One may minimize the free energy functional in the combinatorial 
Thermodynamical Bethe Ansatz approach (TBA) [31 HIE], or one may apply 
algebraic and analytical means for the computation of the partition function 
from the quantum transfer matrix (QTM) [HIE]. 

The TBA approach is based on the string hypothesis and yields an infinite 
set of non-linear integral equations (NLIE). However, it is impractical to solve 
the TBA equations numerically due to the infinite number of equations and 
unknowns. Therefore approximations are required in this approach. 

By means of the quantum transfer matrix approach, a finite set of NLIE 
can be derived exploiting analyticity properties of the quantum transfer ma- 
trix. These equations have been shown to be successful in the description 
of thermodynamical properties in the complete temperature range for many 
important models, like the Heisenberg model [3 El [9] and its spin-S 1 gen- 
eralization [TO], the t — J model [TT|, the Hubbard model [12] and SU(N) 
invariant models for N < 4 |13j . 
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Nevertheless, the standard construction of the quantum transfer matrix 
assumes models with isomorphic auxiliary and quantum spaces. Here we are 
concerned with extensions to more general models with non-isomorphic aux- 
iliary and quantum spaces. Important examples of such systems are mixed 
spin chains. These mixed chains have been extensively studied for low and 
high temperatures by use of the TBA equations and finite size scaling for 
isotropic chains [HI [15]. The dependence on magnetic fields was studied 
[ToT [T71 [TBI [T9] and more recently, also the anisotropic generalization was 
considered [20J. 

Our aim is to propose a construction of the quantum transfer matrix 
by replacing the standard "rotation" of vertex configurations of Boltzmann 
weights by conjugated representations, i.e. by the normal Boltzmann weight 
shifted by the crossing parameter. Having this in mind, we can tackle the 
more general situation where the auxiliary and quantum spaces are not iso- 
morphic. As an application of this idea, we study the generic (Si, S2) case 
of alternating spin chains at finite temperature. 

The paper is organized as follows. In section (2[ we outline the basic 
ingredients of the quantum transfer matrix approach. In section [3[ we define 
the alternating spin chain and its properties. In section HJ we derive the 
set of non-linear integral equations. In section [5j we present our numerical 
findings for the solution of the NLIE. Section [6] is devoted to the calculation 
of the thermal Drude weight for the case Si = S 2 . Our conclusions are given 
in section [3 
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2 Quantum transfer matrix 

We are interested in the computation of the partition function Z = Tr e _/3W 
in the thermodynamical limit, on the condition that Ti is an integrable local 
Hamiltonian derived from some row-to-row transfer matrix. 

In general, transfer matrices can be constructed as ordered products of 
many different local Boltzmann weights £^(A), where A denotes the spectral 
parameter. These weights can be considered as matrices on the space A, 
usually called auxiliary space, which is related to the degrees of freedom on 
the horizontal lines of a two dimensional vertex model. The matrix elements 
of C^i(X) are operators acting non-trivially on the site i of the quantum space 
Yli=i V-i °f a chain of length L and are related to the degrees of freedom on 
vertical lines. 

The product of Boltzmann weights 



defines the monodromy matrix 7^(A). Here we allowed for non-isomorphic 
spaces Vi. This way, A/u(A) - also called ^-operators- may have different 
representations for the L many quantum spaces Cj\i{X) = C A f l \\). The 
labels for different representations, a, Pi, may take for instance integer values 
a, Pi = 0, . . . , L — 1 and C^- a \\) denotes the isomorphic representation. 
Then the row-to-row transfer matrix is the trace over the auxiliary space of 
the monodromy matrix, 



The transfer matrix constitutes a family of commuting operators [T(A), T(/x)] = 
0, provided there is an invertible i?-matrix acting on the tensor product 
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T A (X) = jC A l(\)Cal-iW ■ ■ ■ 4ai(A), 



(1) 



T(A) = Tr^ [T A {\)\ . 



(2) 



A ® A, such that 

J2«(A - n)C^\X) ® £<**V) = ® ^ A) (A)i?W(A - a*). (3) 

In order to have an associative algebra, the .R-matrix is required to satisfy 
the Yang-Baxter equation 

R$(\)R&\\ + aO^SV) = R^(^)R^(X + fi)R&\\). (4) 

The simplest solution of d2J) occurs when auxiliary and quantum spaces 
Vi are isomorphic implying that C^' a \X) = P^R^ (A) , where P i2 is the 
permutation operator. 

The conserved charges are obtained through the derivatives of the loga- 
rithm of the transfer matrix 

.7<»> = ^ln[T(A)]| fco , (5) 

and the Hamiltonian corresponds to the first derivative, TC = J' ( - 1 \ Therefore, 
we can relate the transfer matrix and the Hamiltonian in the following way 

T(A) = T(0)e xn+o{x2 \ (6) 

where T(0) plays the role of a kind of right multiple-step shift operator [13] 
for a general distribution of /^-operators (X). 

Let us consider that in addition to relation ([3]) the ^-operators satisfy 
the following symmetry properties 

Unitarity: C[ a / ] (X) C[ a /\- X) = Ca,/ 3 (A)Id 1 ® Id 2 , (7) 
Time reversal: ^'^(A)* 1 = (X) h } (8) 

Crossing: £ { ^\X) = ^(X)M 1 £^\-X - p/Wf 1 , (9) 
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where ( a ,p(ty an d S a ,pW are scalar functions and p is the crossing parameter. 
Here Id, and U denote the identity matrix and transposition on the i-th. space, 
Mi = M ® M2 where M is some scalar matrix. 

Now, we can define an adjoint transfer matrix T(A) as follows 

L 



7XA) = nW^(A)Tr4 



C 



(»,Pl), 

al y 



-\-p)C 



(a,p L -i), 
AL-l \ 



1=1 



(10) 



and by using the properties (JHJ|9j) we can rewrite the transfer matrix T(A) 
such that, 



T(A) = Tr^ £<g*>(A) . ..C^\X)C^rW 



{*,p L -\), 



*{a,p L ), 



Here we can see that, due to unitarity (JTj), the logarithmic derivative results 
in the same Hamiltonian TC = Ti. and T(0) corresponds to the left multiple- 
step shift operator, such that T(0)T(0) = Id where Af = Yif=i Ca,ft(0). 
In analogy to (jBj), we can write the transfer matrix T(A) as 



T(A) = T(0)e 



XH+0(X 2 ) 



;i2) 



Using ([6]) and ( |T2l) we can rewrite the partition function Z in terms of 
the transfer matrices T(A) and T(A) by considering the Trotter limit, 



lim Tr 
lim Tr 



' e - 2 l-H)N/2 



(T(-r)T(-r)) 



AT/2 



1 



N' 



(13) 
(14) 



The partition function (fT4"l) can be related to a staggered vertex model 



with alternating rows T and T. In this case we need to know all the eigen- 
values of these two transfer matrices to obtain the partition function in a 
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closed form. This is due to the fact that the eigenvalues of both transfer 
matrices depend on the length of the quantum chain L and in particular on 
the Trotter number N, such that for N — > oo all gaps close. However, we can 
circumvent this problem by rewriting (j!4p in terms of the column-to-column 
transfer matrix describing transfer in chain direction and hence is called the 
quantum transfer matrix 



T?™(x) 



...C^\ix + T-p)C^(ix-r)]. (15) 

Each of these objects has a well defined largest eigenvalue separated by a 
gap from the rest of the spectrum, even in the limit N — > oo. Therefore, only 
the largest eigenvalue is required for the computation of the partition func- 
tion. Here x is the spectral parameter associated with the vertical line ensur- 
ing the existence of a commuting family of matrices, T s ®™(x), T- (x 1 ) = 
0. However, of direct physical relevance is x = for obtaining the partition 
function, 

(16) 



lim Tr 

N->oo 



J\f N / 2 ' 



n^™(o) 

Next, we address the identification of the largest eigenvalue of the product 
of the quantum transfer matrices (x). In general, the determination of 

the largest eigenvalue of the product of matrices Yli=i T®™ (x) would require 
the knowledge of all the eigenvalues of all transfer matrices T^™(x), which 
could turn out to be a more involved problem than the staggered model 
mentioned above. 

Nevertheless, this problem can be overcome under certain conditions. 
For instance, for the case of mixed spin chains all of the transfer matrices 
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commute according to the Yang-Baxter equation and the largest eigenvalues 
of the individual transfer matrices correspond to the same eigenvector. This 
implies that the largest eigenvalue of the product of L different transfer 
matrices is nothing than the product of the largest eigenvalues of the quantum 
transfer matrices. In this work, we will restrict to this specific case. 

Here we are interested in the free energy and its derivatives, so we have to 
consider the logarithm of the partition function in the infinite length limit. 
As the eigenvalues Af™(x) depend only on the Trotter number, we can 
first take the infinite length limit and later the infinite Trotter number limit, 
which reads 



Before closing this section, we would like to mention that the properties 
([319]) are also satisfied by many isomorphic self-crossed models [21]. For the 
SU (N) case with N > 2, the property Q reduces to the standard "rotation" 
of the vertex configuration of the Boltzmann weights. 

3 Alternating spin chains 

In the previous section, we used unitarity, time reversal and crossing prop- 
erties to construct the quantum transfer matrix considering general repre- 
sentations of C%i (A). From now on, we consider (for an even number of 
lattice sites L) the alternation of two different representations of the group 
SU{2) with spin S\ at odd sites and spin S 2 at even sites, i.e. /3 2 i— l — Si 




(17) 



(18) 
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and fai = 5*2 . In order to have a Hamiltonian with local interactions we fix 
a to be identical to the spin Si representation (equivalently we could have 
chosen 5*2). 

The monodromy matrix ([1]) becomes 

Tf - S2) (A) = cf L A \\)jO,f L fl\x) . ..4t 52) (A)4t Sl) (A), (19) 

with the auxiliary space A = C 25l+1 , and (A) resp. Cj^ (X) are the 

/^-operators with spin Si representation in the auxiliary space and S 2 resp. 
Si in the quantum space. 

The above SU(2) invariant /^-operators can be obtained through the fu- 
sion process [22J. Its explicit form conveniently normalized is given by 

S1+S2 

4f 1,52) (A)= £ fi(X)Pi, (20) 

J=|Si-S 2 | 

wherJJ /,(A) = nfit+i (Q) n1S(A + S 2 - ^ + j) and P ; is the projector 
onto the SU(2)i in the Clebsch-Gordon decomposition SU(2)s 1 ® SU(2)s 2 - 
This operator is represented by 

S1+S2 o _ o 

p,= n pi) 

k=\S x -S i \ 
k^l 

with x z = I [Z(Z + 1) - Si (Si + 1) - S 2 (S 2 + 1)] and the SU{2) generators 
S a = (S x a ,Sv,S z a ) for a = 1,2. 

The operator fTSOj) is a solution of ([3]) with the following .R-matrix 

i2Sf l) (A) = Pi 2 £g lA) (A). (22) 



lr rhe symbol * shall remind that the possibility j — Si — S2 is excluded throughout 
this work. 
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It satisfies the properties ([7}|H]) with scalar functions given by (si,s 2 W = 
Ttj=i{(£>2 ~ Si + j) 2 ~ an d SSi,s 2 {ty = (— l) 251 and crossing parameter 
p = 1. The matrix M is an anti-diagonal matrix whose non-zero elements 



are M,- 



rf- 52) (A) 



L h3 ~ y *-> J «,25i+2-j- 

The Hamiltonian associated to the transfer matrix T(A) = Tr^ 
has terms with two and three site interactions. Its generic expression is given 
by 



odd i 



-1 



<9A 

L -i-2,i 



(0) 



A=0 

i d 



A=0 



(23) 



where periodic boundary conditions are assumed. For illustration, the Hamil- 
tonian for case S\ — 1/2, S% = S is given explicitly by [15] 

2 

^ ( <?i-i • Si + Si - a i+ i + { ■ Si, Si ■ a i+ i 



2 U+i 



+ ( 5(5+1)) ^5U-<x i+1 



L 

4 



(24) 



One of the consequences of the alternation of two different spins is that 
we have two quantum transfer matrices to work with. We denote them by 
T^ix) and T( s *\x), such as 



J>(Sa)| 



X 



t qtm {x) = TlVa [4^, (ia . + T _ p) 4^) (Lr _ r ) 

. . . 4y + r - p)£gr 5l) (i^ - r)\, (25) 



where the vertical spaces are V a = <C 2Sa+1 and a = 1, 2. 

The transfer matrices (125]) for a = 1, 2 commute due to the Yang-Baxter 
relation [23J . Therefore, they can be diagonalized simultaneously. It can also 
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be deduced from [23] that their largest eigenvalues correspond to the same 
eigenstate. Hence the largest eigenvalue of the product T^ Sl \x)T^ S2 \x) is 
the product of the largest eigenvalues of T^ Sl \x) and T^ S2 \x). 

For the analysis of the spectra we use the fusion hierarchy for the quantum 
transfer matrix T^' (x) , in analogy to the fusion of ^-operators. The algebraic 
relations read (see e.g. [TO] ) 

TV\x)T®(x + i(j + -)) = aj (x)T^\x + -)+ a j+l (x)T^\x - -), 

2 2 2 

T(°\x) = a (x)ld, j = 1,1,1,... (26) 

where dj(x) = nf=i 4>+( x + Ki — Si + I — l))(f)-(x + i(j ~ $i + I)) and 
<f> ± {x) = {x±ir) N ' 2 . 

From the fusion hierarchy with bilinear and linear expressions in T (j26|) . 
one can obtain another set of functional relations [24|, usually called T- 
system, with exclusively bilinear expressions 

T®(x + l -)T^\x - ~) = T u ^\x)T^\x) + fj(x) Id, (27) 

where = Ui=i M% ~ Kj ~Si + l + \))<P-(x - i(j -S 1 + l- \))<j> + {x + 
i(j — S± + I — |))0_(a; + i(j — Si + / + |)) for any j integer or semi-integer. 

Equally important is a set of functional relations referred to as the Y- 
system, which is a consequence of (l27j) . It is written as 

+ l) y 0')(a; - I) = y(^5)( x )y0'+5)( x ) ; (28) 

where yCi>(x) = Z^Ig^W and y(%) = 1 + j/C^x). 

Lastly, we introduce a Zeeman term Ti, = 7i — hS z . This term represents 
the coupling of the magnetic field h to the spin S z = Sii+E i=1 

odd i ' even i ' 
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It can be introduced inside the trace of the partition function such as, 

Z = & Tr [(r(-r)T(-r))" /2 e^] (29) 

Alternatively, it can be considered as a diagonal boundary term on the ver- 
tical lines along a horizontal seam. This redefines only trivially the quantum 
transfer matrix 



T^\x) = TT Va [GaC { " a N Sl \ix + T-p)C iS - Sl) ^ 



VaN-1 



IX 



. . . C^\ix + r - p)4y >\ix - r)], (30) 

where Q a is a diagonal matrix whose non-zero elements are (G a )i,i — e /3h ^ Sa+1 ~ z \ 
The eigenvalues A^(x) associated to T^\x) also satisfy the functional 
relations (l26H28p . This is due to the commutativity property among different 
TV'(x). This way, we obtain the eigenvalues at any fusion level in terms 
of the first level eigenvalue through the iteration of the relations ( |26l) and 
( f27l) . Alternatively, we can proceed along the same lines as [25] applying the 
algebraic Bethe ansatz to the case of twisted boundary conditions. 

In both cases we end up with the eigenvalues of the quantum transfer 
matrix (|3"U1) . 



2j+l 

£ 

m=l 



A '»W = VAf)(x), (31) 



A m [X ) t + , m[ x)t^ ri[ x+i) Q{x _ + . _ m))Q{x _ +j _ m)) 

(32) 

where t<g n {x) = f[ ^ - % ~ f[ 0±(s-i(j-Si+O) and Q(x) = 

i =j - m +2 ( p±y x - 1 y L + ^)) l=1 
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niLi(x — xi). The corresponding Bethe ansatz equations can be written as 

0h <f>+(xi - + \))<t>-(xi - iCgi ~ \)) _ tt xi - Xj - i 
6 0_(^ + i(S 1 + i))0 + (x, + i(S 1 -i)) \\x l -x j ■ V [ > 

According to the previous section, we only need to know the largest eigen- 
value in the limit iV —>■ oo to describe the thermodynamics of the one di- 
mensional quantum model. Then for instance by numerical analysis of the 
Bethe ansatz equation fl33l) for small N we see that the largest eigenvalue 
lies in the sector n = SiN. However, the limit N oo cannot be considered 
numerically. So, we need to encode the Bethe ansatz roots in such a way 
that the free energy can be evaluated independently of the exact knowledge 
of the individual roots. 

One possible way is to define a set of suitable auxiliary functions de- 
pending on the Bethe ansatz roots. Then by exploiting the above and fur- 
ther functional relations we eliminate the explicit dependence on the roots. 
Therefore the Bethe ansatz roots for finite N (including the limit N — > oo) 
become encoded in a finite set of auxiliary functions satisfying certain non- 
linear integral equations. 

Such an analysis was already done for many cases, for instance for the 
spin-1/2 Heisenberg chain 0[8|,[9] and its higher spin extensions [10]. In the 
latter case, the auxiliary functions were taken as a subset of the ^/-functions 
complemented by two "novel" functions which reduce the infinitely many 
functional relations (|28p to finitely many. This is the starting point of the 
next section. 
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4 Non-linear integral equations 



In this section, we introduce a suitable set of auxiliary functions and explore 
its analyticity properties to obtain a finite set of non-linear integral equations. 
These auxiliary functions turn out to describe the largest eigenvalue of (1311 
and consequently the free energy ffl8|) at finite temperature. Specifically, we 
need to define 2s + 1 auxiliary functions, where s = max(Si, S2). We will 
proceed along the lines of [TU] and take as the first 2s — 1 auxiliary functions 
the ^-functions 

a\, x A [j -^(x)A^)(x) 1 1 

y w= — wr^' 3 = 2 s -2- < 34 > 

The two remaining functions are defined as 

«,) = + *) + •■■ + *&>(* + a (35) 

^2s+l\ X + 2) 

- = A^-jH^A^-j) 
*<•>(*- i) 

In addition to this, we introduce a shorthand notation for simply related 

functions B(x) := 1 + b(x), B(x) := 1 + b(x) and Y^\x) := 1 + y^\x) for 

7 = i s - i 
j 2 ' ■ ■ ■ ' 2 ' 

In conformity with the previous definition, we note that B[x) = -^y- 



^7+1(^+5) 

and B[x) = f- with product B{x)B{x) = Y (s >(x). This implies for the 

first (2s — 1) functional relations ff28l) 

y 0) (a; + _ 1) = f or j = i, 1, . . . , s - 1, (37) 

y ('-s)(x + l -)y {s -^{x ~\) = Y {s - 1 \x)B{x)B{x). (38) 
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b(x 



b(x) 



We can write b(x),b(x), B(x) and B(x) explicitly using (1521) such that 

Q(x + i(a + 1)) e /3M^+|)A( s -|) ( z ) 

~ Q(x - is) H?2 + i(s - Si + Z - + i(s - Si + Z + I)) ' 

(39) 

Q(x-i(s + l)) e -^* + ^A(-3)(a;) 



Q(x + is) []S <M* - i(s -Si + l+ \))4>-{x -i{s-S x + l- §)) ' 

(40) 



Q(z - is) 0+(x + i(s - S x + I - \))<f>-{x + i(s-S 1 + l+ §)) ' 

(41) 

Q(x-is) e-"HW(x-i) 



Q(x + is) 0+(x - i(s - Si + Z + - i(s - Si + Z - §)) ' 

(42) 

In this way, it is evident that b(x), b(x) are related to A^ s ~^(x). 

Moreover, A^ s_ 5)(x) is related to Y( s ~^\x) through the definition of y- 
function. This relation can be written as 

A^)( x + i)A( s -D(x - -) = f s _ k (x)Y^\x). (43) 
2 2 2 

At this point, we have a common set of functions which still depend on 
the Bethe ansatz roots and whose limit N — > oo is still to be performed. 
However, this dependence as well as the limit can be worked out easily in 
Fourier space. 

In order to calculate the Fourier transform, we exploit the analyticity 
properties of the eigenvalue of the quantum transfer matrix and the auxil- 
iary functions. Furthermore, these functions should be non-zero and have 
constant asymptotics in a strip around the real axis. This allows us to apply 
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the Fourier transform to the logarithmic derivative of the auxiliary functions, 

?w= < 44 > 

In the cases k < and k > 0, we have chosen a closed contour above 
and below the real axis, respectively. For this reason, it is of fundamental 
importance to analyze the structure of the zeros of the auxiliary functions. 

In particular, the zeros and poles of the auxiliary functions ( |34ll36i) origi- 
nate from the zeros of Q(x) and A^'(x) for j = ~, . . . , s besides those of the 
4>±{x) functions. Therefore, we have to analyze the qualitative distribution 
of the Bethe ansatz roots as well as the zeros of the eigenvalue functions 
A0)(a;). 

It is well known that Bethe ansatz roots form 25'i-strings in the particle 
sector n = S\N. These roots have imaginary parts placed approximately 
at (Si + | — I) for I — 1, . . . , 25*i [23J. Concerning the zeros of A^(x) for 
j = ~, . . . , s, we have verified numerically that their imaginary parts are 
placed at ±(j — Si + I) for / = 1, . . . , 2S\ and I ^ S\ — j. 

By direct inspection of (134113 7H43]) . we note that almost all auxiliary func- 
tions are free of zeros and poles in a strip containing —1/2 < Q(x) < 1/2. 
The exceptions are ^(x) for Si < S2 and b(x),b(x) for Si > S2, which 
should be treated separately. 

This way, the position of the zeros and poles of the auxiliary functions 
depend on the relative magnitude of Si and S2. So, we have to split our 
analysis in three parts: Si < S2, Si = S2 and Si > S2. 



17 



4.1 Si < S 2 

In this case, we have s = S2 in the previous definition. In order to deal with 
the problem involving the function y( Sl \x), we define a related function for 
which the problematic zeros and poles at x — ±i/2 are cancelled, 

+ 2)<M£ - 2) 

Consequently, the 2Si-th equation in ( |37l) becomes 
y< s "(* + !)»<*>(* - 1) = *^Z^±|yC*-i) (l) y(ft + i) W , (46) 



and the functions y( Sl >(x± |) can be transformed as usual according to 
On the other hand, we can apply the Fourier transform to the equation (T45|) . 
once it does not have zeros and poles on the real axis. Thus we are able to 
establish a relation between y( Sl ^ and y( Sl ^ in Fourier space, 

y (Sl) {k) =iNsmh[k(3/N}e^ kl/2 + y (Sl \k). (47) 

Now, applying (jUJ) to the functional relations (I37H43I) and (fltjj) we ob- 
tain after a long but straightforward calculation a set of algebraic relations 
in Fourier space. These relation are given in terms of the transformed aux- 
iliary functions y^(k), b(k), b(k), Y^\k), B(k), B(k) and the unknowns 
A^~^(k), A^(k) and Q(k). We can eliminate the unknowns after some 
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algebraic manipulation. Finally, using (|47p we obtain 









' > 




' Y^)(k) 




y (s ^(k) 




d(k) 




Y^(k) 










+ /C(A;) 















Y^-\)(k) 




b{k) 









B(k) 




V Kk) J 




I J 




K B(k) 


where the kernel JC(k) is a (2S2 + 1) X 


(2S 2 + 1) 


matrix given 


/ 


K(k) 














(48) 



JC{k) 



\ 



K(k) 










k(k) 







K(k) 






K(k) 







k(k) 








k(k) 



K(k) F(k) 






k{k) 

-e- k F(k) 



K(k) -e k F(k) F(k) J 

(49) 



with K(k) 



2 cosh [k/2] ' 2 cosh [k/2] 

As the Trotter number N appears only in d(k), we can take the limit 
N — > 00 straightforwardly, 



-\k\/2 



and d(k) 



2 cosh [fc/2] • 



d(Jfc) 



lim iVsinh [k/3/N] 



ikp 



(50) 



2 cosh [k/2] n^oo 2 cosh [fc/2] ' 

The inverse Fourier transform has been applied to fj48l) followed by an 
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integration over x, resulting in 
^ In 2/5) 





( ^ 




^ lny(5)( x ) ^ 




-(3d(x) 




\nY^\x) 






+ /C* 



















\nB(x) 




v -a ) 




{ In 5 (*) j 



Iny^tx) 

(51) 

In 

where <i(:r) = ^^^] and the symbol * denotes the convolution / * g(x) = 
J^oofix ~ y)9(y)dy- The integration constants ±(3h/2 were determined in 
the asymptotic limit \x\ — > 00. 

The kernel matrix is given explicitly by 



lC(x) 



\ 






K(x) 



















K(x) 





K(x) 















K(x) 





























K(x) 



















K(x) 





K(x) 


K(x) 
















K(x) 


F(x) 


-F(x + i 
















K(x) 


-F(x - i) 


F(x) 



) 

(52) 



where K(x) = and F(x) = \^/2] dk - 

Now, we have to derive an expression for the eigenvalue A.^ S2 \x) in terms 
of the auxiliary functions. It is convenient to define a new function 

k^\x) 



A (S 2 ), 



X 



n?=5 - -5-1+ i))<i>-(x + i(5 2 - 5 X + o) ' 
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(53) 



which has constant asymptotics. For x = and finite N, we have In A^ 2 ^(0) = 

_ AT 

+ j In [A/"^ 2 ] , where we have used the 



1 - 



hiA (52) (0) + £ 2 2ln 
fact that AT = ]1?2( 5 2 -Si + l) 



P i 

52-51 +1 N 
L 



Using the Fourier transformed version of (I41ti42|[53|) . we obtain 



,o \ -\k\{S 2 Sx-\) 25i 

A (2) (k) = ik(3 \ — rL/o1 e ' m + K ( k ) B{k) + B(k) 



2 cosh [Jfe/2] 



i=i 



(54) 



Proceeding as before, we apply the inverse Fourier transform followed by 
an integration over x and the determination of the integration constant. In 
this way, we obtain 



lnA (S2) (a;) = (3e {S2 > Sl) {x) + (K * In BB) (x), 



(55) 



where e^ S2 ' Sl \x) is given by 



6 {S2 ' Sl) (x) 



j2 1 \ ^\w. 

i=i 



2 cosh [jfe/2] 



(56) 



At the point x = 0, we can rewrite this integral in terms of the Euler psi 
function, 



r (S 2 ,Si) 



(o) = V 



S 2 + Si + 1 



if, 



S 2 - Si + 1 



(57) 



2 / V 2 
The contribution of the quantum transfer matrix T^ S2 \0) ( 1301) to the free 
energy is given by (TI8]) 



f^ 2 ^ = lim lnA^(0) + i lim iln^/ 2 ], (58) 

2p N-*oo p N,L— >oo L 



1 1 251 1 

lim lnA( S2 >(0) + - V 



2 S 2 - Si + Z 



(59) 
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Therefore, we can write f( S2 > Sl ") explicitly as 

1 



/ 



(5 2 ,Si) 



2Si 

E 

.1=1 



S 2 -S t + l 



S2 + Si±l\ + fS 2 -St + l 



-—(K*lnBB) (0). 



(60) 



4.2 Si = S 2 



In this case, we note that b(x) and b(x) have zeros at x = ±i/2 which are 
presenting some subtleties. These zeros originate from the factor A^ 1- ^) and 
in principle do not present any problems for the computation of the Fourier 
transform of the logarithmic derivative of ( I39MU1) . The problem arises in 
the Fourier transform of (l4"3l) . which is required to eliminate the unknown 
function A^ Sl ~^\ 



Hence, we define a new function A^ Sl a)( 



A (s i-^(x) 



X ) = + (x-i/2)0-f+i/2) ' WhiCh 

does not have any zeros at x — ±i/2. We apply (j44j) to the functional 
relations (I37H43I) with A instead of A. Then we eliminate the unknowns 
A^ Sl ~^(k) and Q(k) and finally we obtain 



/ f,^)(h\ \ ( n \ 



b{k) 






d(k) 

\ d(k) J 



+ £{k) 



( Y^\k) \ 

f^-k)(k) 
B(k) 
V B(k) J 



(61) 



where the kernel K.(k) with the same structure as (j4"51) . is a (25*1 + 1) x (2S\ + 1) 
matrix. 

Applying the inverse Fourier transform to fl6TT) followed by an integration 
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over x, results in 



/ 



\ 



lny( S l~2)(x) 

In b(x) 
\y \nb(x) J 



( 



\ 



kiB(x) 
y ln.B(x) j 



(62) 



-Pd(x) + $\ 
\ -/3d(x) - (3\ J 
where the (2Si + 1) x (2Si + 1) kernel matrix is given by ( J52l) . 

Finally, the largest eigenvalue A^^O) of the quantum transfer matrix 
can be written in terms of the auxiliary functions in analogy to 
the previous case. We just have to set S2 = Si in all expressions (I53H60P and 
obtain, 



lnA (Sl) (0) = f3 



+ (K*\nBB) (0). 



(63) 



Its contribution to the free energy is given by 

'25i 



c(S u Si) 



^1 . f2S 1 + l\ .(: 



2(3 



(K * In BB) (0). (64) 



4.3 5i > S 2 



For this case, the auxiliary functions as well as the set of non-linear integral 
equations are exactly the same as in the previous case S2 = Si. The only 
difference consists in the way how the largest eigenvalue A^ 2 ^(0) is expressed 
in terms of the auxiliary functions. 

According to the definition of the ^-function, we have an equation similar 
to (|43|) which relates A^ 2 ^(x) and Y^ Sa \x). This relation can be written 
explicitly as 

A (5a) (^ + ^)A (5a) (s ~\) = fs 2 (x)Y {S2 \x). (65) 
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\k\(s 2 -Si+i-±) _ e -\k\(Si-s 3 -i+±) _ e -|fe|(i+o) 



Applying (JUj) to (I65II53I) . we obtain 

25i 25i 

7(ife)= ^ e -|fc|(s a -Si+J-i) _ ^ 

/>(Si-5 2 )+a i<(Si-S 2 )-a 

(67) 

where a = when Si — S2 is an integer number and a = 1/2 when Si — S2 
is a half-integer number. Here, we recall that the possibility I = Si — S2 was 
already excluded in the definition of the £-operator (1201) . 

After performing the inverse Fourier transform and integration over x, we 
obtain 

In A {S2) (x) = (3e {Sl ' S2 \x) + (K * lnF (S2) ) (x), (68) 

with e^ Sl ' S2 \x) = Icosh [k/2] ■ ^ ne P ar ticular point x = 0, e^^'^^x) is 
given by 

^ (0) = ,p±*±i)-,pz*±i). (69) 

Lastly, the contribution to the free energy is written in terms of the 
auxiliary function 

"2Si ' 

S 2 - St + I 



f(S»,Si) = I 
17 2 



— (Jjr*lnr<*>)(0). (70) 

It is interesting to compare e^ Sl,S2 '(x) ( 1691) with the previous cases (I57|63p . 
These expressions can be naturally written in a unified form as follows 

(71) 
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According to (1181) . the free energy of alternating spin chains is described by 
the sum of In A (S2) (0) and In A (Sl) (0). As a result of that, the sum of e^ 1 ' 3 ^ 
and e^ Sl ' Sl ^ is the ground state energy of the quantum Hamiltonian TC^ Sl ' S2 \ 



eo = tp 



S 1 + S 2 + 1 



\Sx - S 2 \ + 1 



[2S 1 + 1 



V 



(72) 



which is in agreement with the results based on the 2S , -string hypothesis for 
the cases S\ = 1/2, S 2 = S [15] and S 2 = Si = S [23J. 

The total free energy is the sum of two pieces / = f( S2 ' Sl ^ + f( s ^< s ^). As 
we have seen, the term f( s ^' s ^ at finite temperature can be written as 



./ 



(S 2l Si) 



Jo 



(S 2 ,Si) 



(K * In BB) (0), ifSt<S 2 

(K * In BB) (0), HSi = S 2 ( 73 ) 

(K*\nY^> s ^) (0), if5i>5 2 , 



where /q 



(S 2 ,Si) 



. Here we have to remind that all 



S2-S1+1 

auxiliary functions, including B(x) and B(x), are different for different cases 
Si < S 2 and Si > S 2 . 

We like to mention that results of an analysis similar to that above were 
published in [26] for the study of single Kondo impurities. In the present 
study of bulk properties of lattice models, the integral equations share some 
algebraic structures with those in [26], but have rather different analytic 
properties with respect to the driving terms. 



5 Numerical results 



In this section, we present the numerical results obtained for the specific heat 
and magnetic susceptibility for the cases S± < S 2 , Si — S 2 and S%> S 2 . 
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(S,S) case 




Figure 1: Specific heat c(T) and x(T) magnetic susceptibility versus temper- 
ature T for S = 1/2,1,3/2,2. 



We have solved numerically the non-linear integral equations by iteration. 
The convolutions have been calculated in Fourier space using the Fast Fourier 
Transform algorithm (FFT). Eventually, we have obtained the free energy as 
a function of temperature and magnetic field. 

Instead of performing numerical differentiations to obtain the derivatives 
of the free energy with respect to temperature and magnetic field, we have 
used associated integral equations for the derivatives of the auxiliary func- 
tions. These integral equations arise from the differentiation of the set of 
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non-linear equations, e.g. with respect to the temperature T. 



(1/2, S) case 




Figure 2: Specific heat c(T) and x{T) magnetic susceptibility versus temper- 
ature T for Si = 1/2 and S 2 = S = 1/2, 1, 3/2, 2. 



Lastly, we have used the relation among the derivatives of the auxiliary 
functions reading 

dT 

d 2 



\nB(x) 
\nB(x) 



1 + b(x) dT 



\nb(x), 



b(x) 



i + b(x) \^r lnKx) ) + W^ nh{x) 



(74) 
.(75) 



dT 2 v ' l + b(x) 
This way, we obtained for each increment in the order of differentiation a 
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new set of linear integral equations, where the lower order derivatives appear 
just as coefficients. 

In Figures [TH, we show the specific heat and the magnetic susceptibility 
as functions of temperature for the particular cases Si = S 2 = S, Si = 
1/2, S 2 = S and Si = S,S 2 = 1/2 for S = 1/2, 1, 3/2, 2, respectively. 

The system shows antiferromagnetic behaviour for the first case Si = S 2 . 
At low temperature c(T) presents a linear temperature dependence and x{T) 
approaches a finite value. For the case S 2 > Si, we have finite magnetization 
Mf = at zero temperature and vanishing magnetic field (T = 0, h — 

+ ) in agreement with [18J. In the other limit (T = + , h — 0), we have zero 
magnetization. This is compatible with the fact that at low temperature and 
zero magnetic field x{T) shows divergent behaviour. For finite (even small) 
magnetic field the system becomes polarized presenting finite magnetization 
associated with a drop of x(T). In the last case, Si > S 2 , the system behaves 
as an antiferromagnet. It has zero magnetization in both limits (T — 0, h — 
+ ) and (T = + , h — 0) in accordance with [T7] . 

For the cases S 2 > Si and S 2 < Si, the models present residual entropy. 
The specific values for this quantity can be extracted from low temperature 



asymptotic solutions of the non-linear integral equations. The results are 

given by S res = \ In [2(S 2 - Si) + 1] and S res = \ In 2S j +2 
z z sm 257+2 

and S2 < Si respectively. The latter case was considered in [17] for (Si 



for 5*2 > Si 



1,S 2 = 1/2) using the TBA approach. There, however, the exact value of 
the residual entropy was left open due to limitations of their method. 
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(S.1/2) case 
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Figure 3: Specific heat c(T) and x(T) magnetic susceptibility versus temper- 
ature T for Si = S = 1/2, 1, 3/2, 2 and S 2 = 1/2 . 



6 Thermal current 

In this section, we are interested in the thermal Drude weight D t h(T) at finite 
temperature. We restrict ourselves to the case Si — S%, where the thermal 
current is related to the second conserved charge (jHJ) of the transfer matrix. 

Specifically, we consider the local conservation of energy in terms of a 
continuity equation. This relates the time derivative of the local Hamiltonian 
ifjj+i to the divergence of the thermal current j , H = —Vj E . Here, the 
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local term H ii+ i stands for 



Hu+i — P{ 



d_ 



C 



(Si,Si) 



(A) 



A=0 



ii+l- 



(76) 



8=1 



As the time derivative leads to the commutator with the Hamiltonian, we 



obtain 



H i>i+1 = i [H, H hl+1 (t)} = -i (jf +l {t) - jf{t)) 
where the local energy current jf is given by 



i [Hi-Ui Hu+i] 



(77) 



(78) 



and the total thermal current is Je = J2i=ijF- 

On the other hand, just by comparing the expression for Je and the 
second logarithmic derivative of the transfer matrix J^ 2 \ we obtain 

L d 2 



Je = J {2) + i 



2dX- 



:C5 1>Sl (A) 



A=0 



(79) 



The transport coefficients are determined from the Kubo formula [27] in 
terms of the expectation value of the thermal current Je, such that 

D th (T) = f (Jl) . 



(80) 



In order to calculate the expectation value (Je), we introduce a new 
partition function Z as, 



Z = Tr [exp(-(m-\ n J {n) )}. 



(81) 



In this way, we obtain the expectation values of J^ through the logarithmic 
derivative of Z, 

d - 2 



<9A 2 



InZ 



A 2 =0 
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(82) 



where we used the fact that the expectation value of the thermal current in 
thermodynamical equilibrium is zero (J~e) = 0. 

To compute the partition function Z, we consider the procedure developed 
in [29]. We rewrite the partition function Z in terms of the row-to- row 
transfer matrix such that 



Z = \im Tr [exp (T{u x ) . . . T(u N )T(0)- N )] . 



Tr 



exp ^jim J2{\nT(ui) - lnT(0)}j 



(83) 



The numbers U\, . . . , are chosen in such a way that the following relation 
is satisfied, 

N 



d 

km 2j{lnT( U/ ) - lnT(O)} = -(3— lnT(x) 
i=i 



8 

•n-l u 

x=o ' " dx n 



+ A n i n - 1 ^lnT(x) 



x=0 

(84) 

We can proceed analogously to section[2]and introduce a quantum transfer 
matrix associated to the partition function Z. Instead of the staggered vertex 
model with alternation in vertical direction between T(— r) and T(— r), we 
have now iV different terms of the form T(0) _1 T(u^) for I — 1, . . . , N. As 
T(O)- 1 = T(0)/JV, we can write T(O)- 1 = [(25i)!]- 2L T(-p). So, we have 
the alternation of T(—p) and T(ui) which is a special case of the previous 
sections. 

Therefore, we can proceed along the same lines as before which is equiv- 
alent to substitute (f>+(x) — > n^li^K^) an d <P-( X ) ~^ YliLi 0o(aO where 
4>i(x) = x — iui and <f>o(x) = x. 

In this way, the partition function can be written in the thermodynamical 
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limit in terms of the largest eigenvalue, 



lim -InZ = lnA(O), 

L^tOO Li 



which is written as 



x=0 



(K*lnBB) (0), 



(85) 



(86) 



where £ (x) = e^ Sl,Sl \x). 

The auxiliary functions B and B satisfy the following set of non-linear 
integral equations 

/ hxy^(x) \ ( \ ( lnY&(x) \ 



\ny( Sl ~^\x) 
In b(x) 
y lnb(x) J 







(-0 + \n£^r)d(x) 

\ (-P + Xn£^l)d(x) J 



\nB(x) 
y ln.B(x) j 



Therefore, the thermal Drude weight is given by, 

2 



D th (T)=f3 2 (J 



r(2)' 



2 



_d_ 



lnA(O) 



A 2 =0 



(87) 



(88) 



In Figure HI we show the thermal Drude weight as function of the tem- 
perature for Si = S2 = S. It exhibits a linear behaviour at low temperatures 
and is proportional to the central charge c = -Mj-. This is in agreement with 
the spin-1/2 case [2"U] . 

Before closing this section, we would like to mention that in the general 
case (Si, $2) the thermal current does not look like a conserved current. In 
this case, we cannot provide exact results for the Drude weight. Neverthe- 
less, we are able to provide an exact description of the second logarithmic 
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T 

Figure 4: Thermal Drude weight D th (T) as function of temperature for S = 
1/2,1,3/2,2. 

derivative of the transfer matrix. However, the physical interpretation of this 
quantity has eluded us so far. 

7 Conclusion 



In this paper we managed to construct the quantum transfer matrix for the 
case of non-isomorphic auxiliary and quantum spaces of interacting spins. 
We considered explicitly the generic (S 1 ,S 2 ) case of alternating spin chains 
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and obtained a finite set of non-linear integral equations. These equations 
were solved numerically for the cases Si < Sa and Si > S%. In this way, 
we obtained the specific heat and the magnetic susceptibility as functions of 
temperature. For the particular case Si = S 2 , we also provided results for 
the thermal Drude weight at finite temperature. 

The system behaves antiferromagnetically for Si > S2 and presents finite 
magnetization in the remaining case Si < S2. Interestingly, for all Si 7^ S2 
we have residual entropy at zero temperature which we were able to evaluate 
exactly. Recently, systems with finite entropy at T = attracted interest 
regarding efficient cooling procedures [31] . 

We expect that our results may be interesting for the study of generic 
mixed spin chains [30J. Another interesting issue deserving investigation is 
the physical interpretation of the second conserved charge for the generic 
case (Si, 5*2 ) and its implications on transport properties. 
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